Attractive Potential around a Thermionically Emitting Microparticle 
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We present a simulation study of the charging of a dust grain immersed in a plasma, considering 
the effect of thermionic electron emission from the grain. It is shown that the OML theory is no 
longer reliable when electron emission becomes large: screening can no longer be treated within 
the Debye-Huckel approach and an attractive potential well can form, leading to the possibility 
of attractive forces on other grains with the same polarity. We suggest to perform laboratory 
experiments where emitting dust grains could be used to create non-conventional dust crystals or 
macro- molecules. 



Introduction: The study of the charging of objects 
immersed in a plasma is a classic problem of plasma 
physics with many applications ranging from space 
problems to dusty plasmas to probe theory for plasma 
diagnostics. Recently, the link between condensed mat- 
ter or high energy-density plasma and strongly coupled 
dusty plasmas P| has renewed the interest in the process 
of charging and shielding of dust in plasmas. 

The interaction between a plasma and an object is due 
to the plasma particles that hit the object surface and are 
captured. In absence of other processes, the higher mo- 
bility of the electrons leads to more electron captures by 
the object, which tends to charge negatively. However, 
in certain conditions, other processes need to be consid- 
ered. For example, if the object immersed in the plasma 
is sufficiently warm, a significant number of electrons can 
be emitted by thermionic effect, altering the balance be- 
tween electron and ion captures, reducing the negative 
charge on the object, or even reversing its sign. An ex- 
ample of this process is given by small objects entering 
the Earth's atmosphere (meteoroids). Recent work has 
shown that the heating of the meteoroids, due to their 
interaction with the atmosphere, can produce a consid- 
erable thermionic emission which can lead to positively 
charged meteoroids 0. 

In the present work, we will consider how electron 
emission changes the process of charging of an object im- 
mersed in a plasma, considering self-consistently charge 
collection on the object and the screening by the sur- 
rounding plasma. We focus particularly on the process 
of thermionic emission, but the results also apply to the 
similar cases of photoemission and secondary emission 
that create a current of electrons emitted by the object. 
Two primary conclusions are reached. First, the process 
of electron emission by the object reduces the charges as 
expected by the orbit motion limited (OML) theory 0. 
However, the quantitative effect of the thermionic emis- 
sion predicted by the OML theory is accurate only for 
small objects. We ffiid that, for objects larger than the 
Debye length, the OML becomes grossly inaccurate. Sec- 
ond, in the presence of thermionic emission when the 
object is charged positively, the screening potential de- 
velops an attractive well. In contrast with the typical 



monotonic behavior predicted by the Debyc-Huckel the- 
ory, we observe a potential well due to the presence of 
an excess of electron charge trapped around the emitting 
object. 

Looking at the literature, we have found experimen- 
tal evidence that heated emissive probes can determine 
a non- monotonic behaviour of the plasma potential 
This is commonly called the virtual cathode, namely a 
region of zero electric field associated with a local excess 
of negative charges. However, as far as we know, the im- 
portance of this mechanism on the charging process oc- 
curring in a dusty plasma and its implications have not 
yet been recognized. The consequences of this behaviour 
of the shielding potential can be considerable since poten- 
tial wells can provide regions of attraction for other ob- 
jects with the same sign of charge. Although the present 
mechanism is not the only instance when particles of the 
same charge immersed in a plasma can attract each other 
(see 0, 01 and references therein), the mechanism pre- 
sented here can be tested experimentally. For example 
photoemission could be used instead of thermionic emis- 
sion in existing experiments. Such experiments would be 
best conducted in microgravity (e.g. on the Alpha space 
station) where other attractive mechanisms (e.g. wake 
field 1^ 1^, ion flow alignment 0) may be less impor- 
tant. 

Charging in Presence of Thermionic Emission: 

We consider a spherical, isolated dust grain of radius a 
immersed in a neutral, unmagnetized plasma consisting 
of electrons and singly charged ions. The grain is station- 
ary, located at r = 0. Ion and electron collisions with the 
neutral gas background are neglected. Electrons and ions 
have different masses, rUf, and m^, and temperatures, Te 
and Ti, respectively. The grain has a surface tempera- 
ture Td and can emit thermionic electrons [W being its 
work function). The characteristic lengths of the sys- 
tem are the electron Debye length Adc, the ion Debye 
length Ad! and the linearized Debye length Aduii , defined 
asl/A2,,„ = l/AL + l/A|,i. 

In the simplest model, neglecting any emissions from 
the dust particle surface, the grain is charged by the sur- 
rounding plasma (primary charging). Initially, electrons 
are more mobile than ions and charge the grain nega- 
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tively by hitting its surface. Hence, the grain acquires a 
negative potential and creates an electric field that repels 
electrons and attracts ions. A dynamical equilibrium is 
eventually reached when the electron current to the dust 
is equal to the ion current. The OML theory provides 
a description of the mechanism, and gives the floating 
potential on the dust, 0(a), as a function of the plasma 
properties. Once (j){a) is known, the electric charge on the 
dust is determined by Qoml = 47reoa (1 + a/Xuvm) 4'[o) 
if one considers a Debye-Huckel potential around the dust 
with screening length given by Aniin- Indeed, the OML 
theory is a good approximation for thick sheaths where 
a <^C \d, but breaks down for a ^ Ad 

The presence of electron emission from the dust (either 
photoelectric or thermionic) affects crucially the poten- 
tial distribution around the dust. In this letter we focus 
on thermionic emission. The starting point for a theoret- 
ical analysis of the thermionic current is the Sommcrfeld 
model of a metal where the energy states are uniformly 
distributed and the free electrons have a Fermi distri- 
bution of probability to occupy a certain energy state. 
We have to distinguish between positively and negatively 
charged dust grains. In fact, when the dust grain is neg- 
atively charged, any electron with energy l/2meV^ > ip 
(ip being the minumum energy required to overcome the 
surface barrier) will be emitted, leading to the following 
thermionic current (IQj 
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known as the Richardson-Dushman expression. When 
the grain is positively charged, the situation is slightly 
different as the electrons have to overcome the floating 
potential as well as the surface barrier. The thermionic 
current for a positively charged grain is [lfl| |: 



IdTT^a^emePT^ / e0(a) 



exp 



W + e0(a) 
kf,. 



(2) 

When thermionic emission is added to the OML frame- 
work, the equilibrium floating potential is established by 
balancing the ion and electron currents from the sur- 
rounding plasma with the thermionic current emitted by 
the dust. 

Simulation Method: To study the charging of a 
thcrmionically emitting dust particle, we have developed 
a PIC code 111)] for a spherical plasma with the station- 
ary grain at the center and an outer radius R. The 
problem under investigation requires special boundary 
conditions. At the outer boundary some particles leave 
the system while others must be injected to represent an 
infinite plasma medium outside the simulation domain. 
The algorithm used to inject the particles is the same 
widely used in the literature 0. At the inner bound- 
ary, the plasma particles reaching the grain surface are 
removed from the simulation and their charge is accu- 
mulated to the central dust grain, affecting its floating 
potential. The same injection method used for the outer 



boundary can be applied also to the thermionic emission 
at the inner boundary, but using the dust temperature 
and not the plasma electron temperature to evaluate the 
distribution function of the emitted electrons. In PIC 
simulations, the emitted electrons are followed accurately 
and the electrons that cannot overcome the dust attrac- 
tion return to the dust; it follows that the emission cur- 
rent must always be computed with Eq. (Q, to avoid 
counting the retarding potential twice. 

We have chosen the parameters of the system accord- 
ing to typical experimental conditions. In particular, we 
consider a Maxwellian plasma with electron temperature 
Te = 1 eV, ion temperature T, = 0.2 eV (TjTi = 5) and 
an outer radius of the system R = 500 /xm. Moreover, 
the plasma far away from the dust grain is Maxwellian at 
rest, with density riao = 6- 10^^ part/m'^. These parame- 
ters correspond to the electron Debye length Ado = 96.0 
/im, the ion Debye length Aoi = 42.9 fim and the elec- 
tron plasma frequency = 4.37T0^s~^. The linearized 
Debye length is Xunn = 39.1 /im. The electron mass is 
chosen with its physical value, but the ion mass is only 
100 times larger. This unphysical choice is common in 
the literature and is required to keep the cost of the sim- 
ulation manageable. All the simulations are made with 
an initial number of particles — Ni — 200000, lo- 



cated on a uniform computational grid with Ng 



200 



cells. The time step is chosen to satisfy the CFL condi- 
tion. At — 10~^^ s. The radial position and two velocity 
components (radial and tangential) are stored for each 
particle during the simulation. 

In the simulations, we start from a uniform Maxwellian 
plasma and let the system relax self-consistently until 
the charge on the dust grain and the shielding potential 
around it reach a steady state. At equilibrium, the dust 
charge fluctuates due to collection of plasma particles 
and we consider the mean value deflned as an average 
over a time interval of 70uj~^ s, which is a sufflciently 
large multiple of the dust charging time to provide a fil- 
ter of the high frequency fiuctuations. Hereafter, "time 
average" will imply an average over the last 70uj~^ of the 
simulation, when the steady state is fully reached. 

Note that all the results presented in the paper are 
obtained in absence of collisions between plasma particles 
and neutrals. Thus our results should be reliable when 
the mean free path for collisions, AcoU, is much greater 
than the electron Debye length, Ado- This requirement 
is well met for example for weakly ionized plasmas: a 
glow discharge with pressure p ~ 10 mtorr, degree of 
ionization x ~ 10~^, density of neutrals Ug ^ 10^° m""^, 
plasma density n ^ 10^^ m""^, leads to Aco1i/Adc ~ 65. 

Results: To validate our simulation tool, we consider 
first the charging of a dust particle in absence of any 
emission. We consider a dust of radius a = 10 /im. Since 
a/ADiin — 0.2, the OML theory should be a good approx- 
imation and we expect our code to agree with theoretical 
predictions. This is indeed true, the time average floating 
charge is Qd = — 1.54-10^^^ C, while the one predicted by 
the OML theory is Qqml = -1.64- lO^^^ C. The relative 
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difference is defined as \Qd - QomlI/ max(Qd, Qoml) = 
6%. At dynamical equilibrium, the floating potential is 
(f)d = —1.1386 V, in good agreement with the one given 
by the OML theory, 0oml = -1.1776 V (relative dif- 
ference 3%). Furthermore, the shielding potential fol- 
lows closely the Debye-Huckel expression with a screen- 
ing length equal to the linearized Debye length Aduii- 

We have also considered the primary charging mecha- 
nism for a dust of radius a = 80 /im. Since a/Aoun — 2, 
we expect the OML theory to be unreliable. In fact, 
our code gives Qd ~ —2.04 • 10"^'' C while Qoml = 
—3.19 ■ lO^^"' C with a relative difference of 36%. On the 
other hand, the value of the floating charge defined by 
Qdc = 47reoa (1 -I- o/Ado) <t>{o) is a good estimate of Qd 
{Qoe = -1.92-10"^'' C). This is a consequence of the fact 
that when the dust radius grows, the screening length is 
determined by the electrons 0- The profile of the time 
average shielding potential follows the Dcbyc-Huckel ex- 
pression but now with a screening length equal to the 
electron Debye length Aoe, as predicted in Ref. The 
time average floating potential obtained by the simula- 
tion is 4>d = —1.2044 V and the relative difference with 
respect to (j)oML = -1.1776 V is 2%. Thus, the OML 
theory gives a good estimate also when a/ADiin — 2, pro- 
vided that the screening length is determined by the elec- 
trons. The value of (j)d in the present case is more negative 
than in the case of o = 10 ^m due to the development 
of an absorption barrier that diminishes the ion current 
to the dust. Furthermore, the sheath is wider, of the or- 
der of several linearized Debye lengths. In summary, our 
code has confirmed all the theoretical predictions from 
the OML theory regarding non emitting dust particles. 

Next, we include thermionic emissions. We consider a 
dust at Td = 0.1 eV and with work function W = 2.2 eV 
(representative of some metallic oxides). These param- 
eters lead to a positively charged grain. We have per- 
formed a number of simulations varying the dust radius 
a. Here we focus on the case a = 80 fim to point out the 
most relevant aspects of the role of thermionic emission in 
the charging mechanism. The dynamical equilibrium is 
reached in approximately 2 electron plasma periods, be- 
ing determined essentially by the electron and thermionic 
currents. (The ion dynamics is relevant on a longer time 
scale, creating the non-monotonic behaviour of the ion 
density explained below.) The equilibrium charge is: 
Qd = 8.51 ■ 10"^^ C, where the OML theory predicts 
QoML — 5.17 • 10^^^ C (based on expression (0) and 
ADiin) or Qbc = 3.12 • 10~^^ C (based on expression Q 
and Ado). However, when the thermionic effect is taken 
into account a comparison of the floating charge of the 
simulation and of the OML theory is no longer correct. 
In fact, the numerical factor that dcflncs the floating 
charge from the floating potential depends on the poten- 
tial distribution around the dust which, when thermionic 
emission is present, is not well represented by the Debye- 
Huckel potential (either with Aoun or Adc). Focusing 
on the floating potential, we find (f)d = 0.1016 V while 
(f'OMh = 0.1911 V. Interestingly, when the thermionic ef- 
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FIG. 1: Thermionic emission. Shielding potential (/f>(r) at dif- 
ferent times (left) and time average shielding potential (right) 
as a function of the dust temperature. 



feet is present, the OML theory does not produce an ac- 
curate estimate of the charging mechanism of large grains 
(also in cases when, in absence of thermionic emission, its 
predictions are acceptable). Moreover, we have checked 
that, for small objects, the OML theory is still reliable 
when electron emission is included. 

Figure Q (left panel) shows the shielding potential at 
5 different times of the simulation. A potential well is 
present. The presence of such well is of considerable in- 
terest since it can lead to attractive forces on another 
dust particle, even when it has the same charge. On the 
right panel of Fig. ^ one can see the time average shield- 
ing potential (solid line) and the potential well is clearly 
visible. We have also shown the time average shielding 
potential obtained by two simulations with Td = 0.01 eV 
(dotted line) and Td — 0.2 eV (dash-dotted line). 

Discussion: How can a potential well form around 
the dust? The explanation comes from the comparison 
with the case of absence of thermionic emission. 

In presence of sufficient thermionic emission (which is 
the case considered above) the dust is positively charged. 
The resulting electric field attracts electrons and repels 
ions. As a consequence, the electrons emitted from the 
dust are slowed down in a region very close to the dust. 
The more energetic electrons can escape and contribute 
to the thermionic current emitted from the dust but 
the rest of the thermionic electrons form an electron 
cloud. The electron cloud determines an excess of nega- 
tive charge and leads to an (equilibrium) potential well. 

To support this explanation. Fig. |21 shows a compari- 
son of the time average ion and electron density for the 
case with only primary charging (left two panels) and 
with thermionic emission (right two panels). The densi- 
ties are normalized with respect to the unperturbed ion 
density, pi^oo- For the primary charging case, the densi- 
ties are perturbed roughly to a distance of 4 Adiiii from 
the dust grain; the electron density decreases while the 
ion density increases towards the grain. The last result is 
due to spherical geometry and to the ion angular momen- 
tum j9j. In fact, there are many ions with high angular 
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FIG. 2: Time average normalized electron and ion densities pe 
and Pi for a = 80 ^m: primary charging (left) and thermionic 
emission (right). 



momentum that do not strike the grain, thus leading to 
an increment of ion density in the sheath with respect 
to the equihbrium value. Consider next the case with 
thermionic effect (right panels). The electron density in- 
creases close to the grain both due to thermionic emission 
and to the attractive potential on the dust. On the other 
hand, since the grain is positively charged, the ion den- 
sity diminishes approaching it. It can be noticed that the 
ion density increases from the dust grain somewhat up to 
4 ADiin, reaches a maximum and decreases to the value at 
rest. Clearly, Fig.|3shows the excess of electrons needed 
for the formation of the attractive well observed above. 

We have performed another simulation where we have 
kept the potential on the dust fixed at the same potential 
observed in the simulation described above, in presence of 
thermionic emission. In the present case, we impose the 
dust potential and we are not allowing any thermionic 
emission: only the primary charging is in effect. Note 
that this simulation is actually a description of the well 
known Langmuir probe used in experimental plasma di- 
agnostics. For this case, Fig. |31 shows the potential dis- 
tribution at different times (left panel) and the time av- 
erage shielding potential (right panel). One can see that 
the time average potential around the dust is a decreas- 



ing monotonic function of radius and vanishes asymptot- 
ically. Clearly, this confirms that the excess of electrons 
seen in Fig. El depends on the thermionic electrons. 

Finally, we have found the dust temperature critical 
for the well formation, as shown in Fig. ^ (right). When 
Td = 0.01 eV, the thermionic current is sufficiently small 
so that the grain is negatively charged and the shielding 
potential is an increasing monotonic function of radius 
(as for the primary charging case). When the dust tem- 
perature grows, for example up to Td = 0.1 eV, an elec- 
tron cloud forms around the dust and causes a potential 
well. However, if Td is increased further (for example 
Td = 0.2 eV), more thermionic electrons have enough ki- 
netic energy to escape the dust attraction. At this high 
dust temperature, the local excess of electrons is reduced, 
the potential well disappears and the shielding potential 
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FIG. 3: Langmuir probe. Shielding potential (j}{r) at different 
times (left) and time average shielding potential (right), for 
a dust particle with a = 80 pm. The Debye-Huckel poten- 
tial with Aoiin (dashed line) is shown for comparison. The 
potential on the dust is held fixed at (j}{a) = 0.115 V and no 
thermionic electrons are emitted. 



becomes a monotonically decreasing function of radius. 
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